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AUTOMATIC FREQUENCY CONTROL METHOD AND DEVICE AND DEMODULATOR 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to an automatic frequency 
control method and device, which are applied to digital 
demodulation processing in satellite communication, mobile 
satellite communication, mobile terrestrial communication, 
and a demodulator to which said device is applied. 

2 . Description of the Related Art 

Recently, in satellite communication, mobile satellite 
communication, and mobile terrestrial communication, digital 
modulation and demodulation have been actively researched. In 
particular, in the environment of mobile communications, a 
radio signal is greatly influenced by fading. Therefore, 
various demodulation systems which stably operate even in such 
a fading environment have been considered. Among these systems, 
a system, which is constructed so that absolute coherent 
detection can be executed even in a fading environment by 
estimating and compensating fading distortion by using known 
signals, has been recognized. In the case where fading 
distortion is estimated and compensated after quasi-coherent 
detection or the like is carried out by this system, to estimate 
and compensate fading distortion with high accuracy, it is 
necessary that the frequency offset between the carrier wave 
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frequency of a radio transmission signal and the oscillation 
frequency of a reference signal for quasi-coherent detection 
is small. 

However, in the case where the frequency stability and 
accuracy of the oscillation circuit of a transmitter-receiver 
are insufficient, unless the frequency of a radio receiving 
signal is automatically controlled by eliminating this 
frequency offset by some processing, fading distortion cannot 
be estimated and compensated with high accuracy. 

In mobile communication, transmission and receiving are 
carried out between a fixed station and a mobile station or 
between mobile stations. Therefore, when two stations 
relatively move, the frequency of a radio receiving signal 
deviates due to Doppler variation. Therefore, even if the 
stability and accuracy of the oscillation circuit of the 
transmitter-receiver are good, frequency offset occurs between 
the frequency of the radio receiving signal and the oscillation 
frequency of the reference signal. 

A technique for compensating frequency offset is 
disclosed in, for example, Japanese Patent Laid-Open No. 
93302/1997 titled "DIGITAL MOBILE RADIO COMMUNICATION METHOD" . 
By the conventional technique disclosed in this document, 
frequency offset is eliminated by using phase fluctuation 
information of known signals (pilot signals). 

In this conventional technique, a radio transmission 
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signal into which two-symbol known signals have been inserted 
every inserting period Np from the transmission side is 
transmitted. On the other hand, at the receiving side, the phase 
change amount between the sequential two-symbol known signals 
is calculated, and in accordance with the calculated phase 
change amount, the phase of the radio receiving signal is 
rotated. Frequency offset is thus eliminated from the radio 
receiving signal. 

In electric wave transmission channels, there is a Rician 
fading transmission channel in which direct waves and 
multi-path waves are mixed. In this case, the direct wave is 
Doppler-shifted. Therefore, the frequency fj, of the direct wave 
further deviates by Doppler-shif t amount f^p caused by the 
Doppler-shif ting from the offset f^rsT caused by stability of 
the oscillation circuit as shown in Fig. 24A. 

On the other hand, in the abovementioned prior art, the 
phase change amount between the sequential two-symbol known 
signals is supposed as frequency offset. The frequency offset 
calculated in this case is equivalent to the difference between 
the oscillation frequency f^ at the receiver side and the 
frequency fo of the direct wave. That is, in the prior art, as 
shown in Fig. 24B, frequency control is carried out so that the 
frequency f^ of the direct wave almost coincides with the 
oscillation frequency f^,. In this case, the center frequency 
fn of the Doppler spread deviates from the oscillation frequency 
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fo by the Doppler-shift amount f^p. Therefore, the Doppler spread 
apparently further broadens, so that the frequency 
corresponding to the end of the Doppler spread greatly deviates 
from the oscillation frequency f^. Accordingly, frequency 
offset cannot be satisfactorily compensated. Therefore, bit 
error rate characteristics (hereinafter, referred to as BER 
characteristics) after the radio receiving signal is 
demodulated may deteriorate. 

SUMMARY OF THE INVENTION 

The object of the invention is therefore to solve the 
abovementioned technical problems and provide an automatic 
frequency control method and device by which satisfactory BER 
characteristics can be realized without being influenced by a 
Doppler-shifted direct wave in a Rician fading environment 
where the direct wave has Doppler-shifted. 

Another object of the invention is to provide a 
demodulator by which demodulation accuracy can be improved by 
applying the abovementioned automatic frequency control device 
to this demodulator. 

In order to achieve the abovementioned objects, an 
automatic frequency control method of the invention, which 
controls the frequency of a radio receiving signal by 
compensating frequency offset of the radio receiving signal 
that periodically includes a plurality of adjacent known 
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signals,, comprising, based on the distortion amount of each 
known signal included in the radio receiving signal, the 
frequency of the direct wave of the radio receiving signal and 
the center frequency of the Doppler spread of the radio 
receiving signal are estimated, and based on both of the 
frequencies, frequency offset of the radio receiving signal is 
compensated. 

In this construction, by compensating the frequency 
offset based on the frequency of the direct wave, the 
compensating range of frequency offset can be sufficiently 
secured, and by compensating the frequency offset based on the 
center frequency of the Doppler spread, satisfactory BER 
characteristics can be secured. That is, compatibility with 
securing of a sufficient frequency offset compensating range 
and securing of satisfactory BER characteristics can be 
realized. 

In addition, the abovementioned construction is applied 
to a demodulator, whereby demodulation processing can be 
executed for a radio receiving signal in which frequency offset 
has been satisfactorily eliminated. Therefore, the 
demodulation quality can be improved. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a conception diagram showing the entire 
structure of the radio communication system to which an 
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automatic frequency control method relating to Embodiment 1 of 
the invention is applied; 

Fig. 2 is a drawing showing the format of a burst signal; 

Fig. 3 is a flowchart for explanation of demodulation 
processing; 

Fig. 4 is a drawing for explanation of estimated timing 
of a frequency offset; 

Fig. 5 is a flowchart for explanation of a frequency 
offset estimation processing; 

Fig. 6 is a graph showing BER characteristics; 

Fig. 7 is a flowchart for explanation of the first 
frequency offset estimation processing; 

Fig. 8 is a flowchart for explanation of the second 
frequency offset estimation processing; 

Fig. 9 is a conception diagram for explanation of the 
second frequency offset estimation processing; 

Fig. 10 is a graph for explanation of BER based on 
Embodiment 1 ; 

Fig. 11 is a flowchart for explanation of the first 
frequency offset estimation processing relating to Embodiment 

2 of the invention; 

Fig. 12 is a flowchart for explanation of the second 
frequency offset estimation processing relating to Embodiment 

3 of the invention; 

Fig. 13 is a flowchart for explanation of the second 
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frequency offset estimation processing relating to Embodiment 

4 of the invention; 

Fig. 14 is a flowchart for explanation of the first 
frequency offset estimation processing relating to Embodiment 

5 of the invention; 

Fig. 15 is a flowchart for explanation of the first 
frequency offset estimation processing relating to Embodiment 

6 of the invention; 

Fig. 16 is a flowchart for explanation of frequency offset 
estimation processing relating to Embodiment 7 of the 
invention; 

Fig. 17 is a flowchart for explanation of demodulation 
processing relating to Embodiment 8 of the invention; 

Fig. 18 is a flowchart for explanation of demodulation 
processing relating to Embodiment 9 of the invention; 

Fig. 19 is a flowchart for explanation of demodulation 
processing relating to Embodiment 10 of the invention; 

Fig. 20 is a flowchart for explanation of the first AFC 
processing; 

Fig. 21 is a flowchart for explanation of the second AFC 
processing; 

Fig. 22 is a flowchart for explanation of demodulation 
processing relating to Embodiment 11 of the invention; 

Fig. 23 is a flowchart for explanation of demodulation 
processing relating to Embodiment 12 of the invention; and 
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Fig. 24 is a conception diagram for explanation of the 
relationship of the frequency offset estimating range, 
estimation accuracy, and possible frequency offset. 

DETAILED DESCRIPTION OF THE PREFFERED EMBODIMENTS 

Hereinafter, Embodiments of the invention are explained 
in detail with reference to the attached drawings. 

Embodiment 1 

Fig. 1 is a block diagram showing the construction of a 
radio communication system to which an automatic frequency 
control method relating to Embodiment 1 of the invention is 
applied. This radio communication system comprises transmitter 
1 and receiver 10, and has a function for estimating and 
eliminating frequency offset of a burst signal radio- 
transmitted from the transmitter 1 by receiver 10. 

Concretely, for this radio communication system, a 
satellite communication system, mobile satellite 
communication system, and mobile terrestrial communication 
system can be applied. In the satellite communication system, 
transmitter 1 and receiver 10 are applied to earth stations set 
at different positions on the earth. In the mobile satellite 
communication system, the transmitter 1 is applied to either 
an earth station or a mobile station that is set on the earth, 
and the receiver 10 is applied to an opposite station of the 
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transmitter 1 of the two stations. In the mobile terrestrial 
communication system, the transmitter 1 is applied to either 
a base station or a mobile station, and the receiver 10 is applied 
to an opposite station of the transmitter 1 of the two stations. 

As a mobile station in the mobile satellite communication 
system, a satellite mobile phone of a single mode or a dual mode 
can be applied. In addition, a cellular phone can be applied 
to the mobile station in the mobile terrestrial communication 
system. 

This coiranunication system uses the TDMA (Time Division 
Multiple Access) as a communication method. Therefore, the 
transmitter 1 radio-transmits a burst signal synchronized with 
a predetermined time slot to the receiver 10. The receiver 10 
demodulates the receiving burst signal, eliminates frequency 
offset, and then restores the original signal. 

The transmitter 1 comprises known signal insertion 
circuit 2, modulation circuit 3, and transmission circuit 4. 
The known signal insertion circuit 2 periodically inserts known 
signals into information signals to be transmitted, and 
generates burst signals before modulation. Concretely, as shown 
in Fig. 2A, the known signal insertion circuit 2 inserts 
Np-symbol sequential known signal rows (hereinafter, referred 
to as known signal blocks) into (Nj,-Np) -symbol information at 
every Nj.-symbol periods. Np is an integer of 2 or more {NpS:2). 
The timings at which the known signals are inserted are shown 
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by (kNp+i)Ts. Herein, k indicates the insertion order of the 
known signals, i is a value between zero and (Np-1) (0:£i^Np- 
1). Tg is a symbol period. 

The known signal insertion circuit 2 provides a burst 
signal before being modulated to the modulation circuit 3. The 
modulation circuit 3 modulates this non-modulated burst signal, 
and outputs it as a burst signal after being modulated. The 
modulation circuit 3 synchronizes the modulated burst signal 
with a predetermined time slot and provides the signal to the 
transmission circuit 4. Concretely, as shown in Fig. 2B the 
modulation circuit 3 synchronizes burst signals B, B+1, B+2, 

B+3 with predetermined time slots SI, S2, S3, S4 and 

provides the signals to the transmission circuit 4. The 
transmission circuit 4 superposes the burst signals on electric 
waves and transmits them to the receiver 10. 

The receiver 10 includes receiving circuit 11 and 
demodulator 12. The receiving circuit 11 is comprised of an 
amplifier, a frequency converter, and the like, and converts 
the frequency of the receiving burst signal which is a radio 
receiving signal into an intermediate frequency and outputs it 
as a receiving IF signal. The demodulator 12 demodulates the 
receiving IF signal outputted from the receiving circuit 11 by 
quasi-coherent detection, and restores the original 
information signal. More concretely, the demodulator 12 
comprises frequency converter circuit 21, two A/D 
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(analog/digital) converter circuits 25a and 25b, and DSP 
(Digital Signal Processor) 26 as an automatic frequency control 
device or digital signal processing device, and restores the 
original information signal by digitally executing 
demodulation processing by the DSP 26. 

The frequency converter circuit 21 converts a receiving 
IF signal periodically including a plurality of known signals 
into an analog baseband signal, and the circuit includes one 
oscillation circuit 22 , two multiplication circuits 23a and 23b, 
and jr/2 phase shift circuit 24. The oscillation circuit 22 
generates local oscillation signals with a predetermined 
oscillation frequency. The local oscillation signals generated 
by the oscillation circuit 22 are inputted into two 
multiplication circuits 23a and 23b. In this case, the local 
oscillation signal is inputted into the multiplication circuit 
23a for I channel via the jt/2 phase shift circuit 24. The jt/2 
phase shift circuit 24 shifts the phases of the local 
oscillation signals by 3t/2. Therefore, the local oscillation 
signals whose phases are shifted by n/2 from each other are 
inputted into the multiplication circuits 23a and 23b. 

The multiplication circuits 23a and 23b each mix the 
receiving IF signals and the local oscillation signals. As a 
result, I channel and Q channel analog baseband signals are 
generated. These generated analog baseband signals are supplied 
to the A/D converter circuits 25a and 25b respectively. The A/D 
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converter circuits 25a and 25b each convert the analog baseband 
signals into I channel and Q channel digital baseband signals. 
The digital baseband signals are supplied to DSP 26 
respectively . 

The DSP 2 6 receives the inputs of the digital baseband 
signals generated by the A/D converter circuits 25a and 25b, 
and based on the distortion amounts of the known signals 
included in the inputted digital baseband signals, estimates 
the frequency of the direct wave and the center frequency of 
the Doppler spread, and compensates frequency offset based on 
these estimated frequencies. Thus, the DSP 2 6 automatically 
controls the frequencies of the digital baseband signals. 
Furthermore, the DSP 26 eliminates fading distortion from the 
digital baseband signals from which frequency offset has been 
compensated, and then demodulates the digital baseband signals, 
whereby the original information signals are restored. 

More concretely, the DSP 26 has storage 2 6a comprised of 
a ROM or the like. In the storage 26a, a demodulation processing 
program which is a computer program is stored. By executing the 
demodulation processing program stored in the storage 26a, the 
DSP 26 executes serial demodulation processing including 
estimation of the abovementioned two frequencies, compensation 
for frequency offset, compensation for fading distortion and 
demodulation . 

Fig. 3 is a flowchart for explanation of the demodulation 
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processing to be executed by the DSP 26. The DSP 26 applies 
filtering such as waveform shaping processing to the digital 
baseband signals supplied from the A/D converter circuits 25a 
and 25b (step SI). Thereby, the DSP 26 eliminates the high 
frequency components such as noise components higher than a 
predetermined cutoff frequency from the digital baseband 
signals . 

Thereafter, the DSP 26 executes the Nyquist point 
detection processing (step S2). Concretely, the DSP 26 detects 
the Nyquist point of the digital baseband signal to obtain 
digital baseband signal r(kNp+i) corresponding to the Nyquist 
point. The digital baseband signal r(kNp+i) corresponding to 
the Nyquist point is expressed as Formula (1) below. In Formula 
(1) below, c(kNp+i) shows a distortion amount caused by fading. 
In addition, A shows the amplitudes of the signals, and b(kNp+i) 
shows symbol values. Furthermore, n(kNp+i) shows gauss noises. 
r(kNj. + i) = Ac(kNp + i)b(kNp + i) + n(kNp + i) ( i ) 
Thereafter, the DSP 26 executes automatic frequency 
control processing (Step S3 through S5). In the automatic 
frequency control processing, frequency offset based on the 
oscillation frequency of the oscillation circuit 22 as a 
reference is estimated, and the estimated frequency offset is 
eliminated from the digital baseband signal r(kNp+i), whereby 
the frequency of the digital baseband signal r(kNp+i) is 
automatically controlled. In other words, the DSP 26 estimates 
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the frequency of the direct wave and the center frequency of 
the Doppler spread based on the distortion amounts of the 
plurality of known signals periodically included in the digital 
baseband signal corresponding to the radio receiving signal, 
compensates the frequency offset of the digital baseband signal 
r( kNp+i) based on both the abovementioned frequency, and thereby 
automatically controls the frequencies of the digital baseband 
signal r(kNp+i) . 

Further in other words, the DSP 26 estimates frequency 
offset of the digital baseband signal r(kNp+i) from the 
frequency of the direct wave and the center frequency of the 
Doppler spread that are estimated based on the distortion 
amounts of the plurality of known signals periodically included 
in the digital baseband signal r(kNF+i) corresponding to a radio 
receiving signal, eliminates the estimated frequency offset 
from the digital baseband signal r (kN^+i) , whereby the frequency 
of the digital baseband signal r(kNp+i) is automatically 
controlled. 

More concretely explaining the automatic frequency 
control processing, the DSP 26 executes frequency offset 
estimation processing first (step S3). More specifically, the 
DSP 26 estimates a phase difference es(mLNj.) corresponding to 
the phase rotation amount of one symbol as frequency offset 
based on the digital baseband signal r(kNy+i). Still more 
specifically, the DSP 26 determines a phase difference 6s(mLNp) 



15 



as frequency offset at each estimating time mT (m: natural 
number) that comes after each estimating period T (=LNj,Ts) (see 
Fig. 4 ) based on the distortion amounts of the plurality of known 
signals included in the digital baseband signal r(kNp+i). 

Next, the DSP 26 executes integration processing for 
integrating the determined phase difference 6s(mLNp) by one- 
symbol period Tg (step S4). Concretely, the DSP 26 cyclically 
sums phase differences es(mLNp) at each symbol period Tg as shown 
in Formula (2) below. Thereby, the DSP 26 obtains an accumulated 
phase difference 6 ( (m-1 )LNp+lNj,+i) . In Formula (2) below, 1 
indicates identification numbers corresponding one to one known 
signal blocks, and is a value of 0 or more and (L-1) or less 
(0:slsL-l). L shows the number of known signal blocks included 
in the estimating period T. 

e( (m - IJLNj. + iNj, + i) = e( (m - l^LNj. + iNp + 1 - 1) + eg(nLLNp) 

(2) 

Thereafter, the DSP 26 executes frequency offset 
elimination processing for eliminating frequency offset from 
the digital baseband signal r(kNF+i) based on this accumulated 
phase difference e( (m-1 )LNF+lN5,+i) (step S5). Concretely, the 
DSP 26 rotates in reverse the phase of the digital baseband 
signal r(kNp+i) by an amount in accordance with the accumulated 
phase difference e( (m-1 )LNp+lN5.+i) . Thereby, frequency offset 
can be eliminated from the digital baseband signal r(kNs,+i). 
That is, as shown in Formula (3) below, the DSP 2 6 obtains a 
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digital baseband signal rR(kNp+i) in which frequency offset has 
been eliminated. Thus, the DSP 2 6 automatically controls the 
frequency of the digital baseband signal. 
r^(kN^ + i) = rj,( (m - l)LNp + IN^ + i) 

= r( (m - l)LNr + iNp + i) exp[-je( (m - IJLN^ + IN^ + i) ] 

(3) 

Thereafter, the DSP 26 executes fading distortion 
compensation processing for estimating and eliminating fading 
distortion from the digital baseband signal rR(kN5,+i) (step S6) . 
Concretely, the DSP 26 extracts Np-symbol known signals from 
the digital baseband signal rR(kNj,+i). The frequency offset 
caused by the transmitter 1 and receiver 10 has been already 
eliminated from these extracted known signals. 

The DSP 26 detects fading distortion based on the 
extracted Np-symbol known signals. The DSP 26 executes 
interpolation processing such as gauss interpolation, Wiener 
interpolation, based on the detected fading distortion. The DSP 
26 thus estimates fading distortion in the information signal. 
Furthermore, the DSP 26 eliminates the estimated fading 
distortion. Thereby, the DSP 26 executes fading interpolation 
for the digital baseband signal rR(kNp+i). 

Thereafter, the DSP 26 executes data decision processing 
(step S7). Concretely, the DSP 26 dicides the original 
information signal based on the digital baseband signal in which 
fading has been compensated. Thus, the demodulation processing 
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is completed. 

Fig. 5 is a flowchart for explanation of the frequency 
offset estimation processing. This frequency offset estimation 
processing includes the first frequency offset estimation 
processing and the second frequency offset estimation 
processing. That is, the DSP 26 estimates frequency offset by 
combining the first frequency offset estimation processing and 
the second frequency offset estimation processing. 

In the first frequency offset estimation processing, the 
frequency of the direct wave of the radio receiving signal is 
estimated as the first frequency offset. In the second frequency 
offset estimation processing, the center frequency of the 
Doppler spread of the radio receiving signal is estimated as 
the second frequency offset. The DSP 26 estimates the last 
frequency offset from the first and second frequency offsets 
estimated by these two frequency offset estimation processings . 

Thus, the reason why the two frequency offset estimation 
processings are combined is as follows. As explained in the 
subsection "Description of the Related Art", in the case where 
the frequency f^ of the direct wave is regarded as frequency 
offset and compensated, as shown in Fig. 24B, the Doppler spread 
apparently further broadens, and the BER characteristics 
deteriorate. On the other hand, in the case where the center 
frequency fj^ of the Doppler broadening is regarded as frequency 
offset and compensated, as shown in Fig. 24C, the direct wave 
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has frequency deviation corresponding to the Doppler shift 
amount f^p, however, the Doppler spread is left at the original 
degree of broadening. Therefore, in this case, the 
deterioration of the BER characteristics is less than that in 
the case where the frequency f^ of the direct wave is regarded 
as frequency offset and compensated. 

Fig. 6 is a graph showing the BER characteristics. As can 
be clearly understood from this graph, at the point where the 
compensation for a frequency close to the center frequency f^ 
of the Doppler spread suppresses the deterioration of the BER 
characteristics more greatly than in the case where the 
frequency f^ of the direct wave is compensated, it is understood 
that the center frequency f^i is the optimum frequency offset. 
That is, in terms of prevention of deterioration in the BER 
characteristics, it can be said that the center frequency 
of the Doppler spread being regarded as the frequency offset 
is better than the frequency f^ of the direct wave being regarded 
as the frequency offset. 

On the other hand, in the case where the center frequency 
fji of the Doppler spread is regarded as frequency offset and 
compensated, prevention of deterioration in the BER 
characteristics is possible, however, there is a problem where 
the frequency offset compensating range is relatively narrow. 
On the other hand, in the case where the frequency f^ of the 
direct wave is regarded as frequency offset and compensated. 
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a relatively wide frequency offset compensating range can be 
secured. 

Concretely, on supposition that the phase rotates by AQ 
in a certain time At, the frequency offset Af can be expressed 
by Formula (4) below. 

Af = ( 4 ) 

2jt At ^ ' 

When this frequency offset Af is estimated by using the 
phase rotation amount between the known signals, the offset can 
be expressed by Formula (5) below. In Formula (5) below, AOp 
shows the phase rotation amount between the known signals. In 
addition, Rs(symbol/s) shows the transmission rate of the 
signals . 

1 A0 1 Rg 

Af ^AOp (S) 

The detecting range of the phase rotation amount AOp is 
-jt^Ap6:sji, and as a result, the estimating possible range of the 
frequency offset Af can be expressed by Formula (6) below. 

2Nf 2Nf ^ ' 

When the frequency f^ of the direct wave is estimated, 
since the phase rotation amount between the adjacent known 
signals is used, the f,, becomes equivalent to the value obtained 
by Formula (6) upon substituting 1 for Np. That is, the frequency 
offset estimating range fpETi in the first frequency offset 
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estimation processing becomes a range between -Rs/2 and Rs/2. 

In addition, when the center frequency f^, of the Doppler 
spread is estimated, since the phase rotation amount between 
the known signal blocks inserted for each Np symbol is used, 
the frequency offset estimating range f^ETj becomes the same as 
the range expressed by the abovementioned Formula (6) . That is, 
the frequency offset estimating range f^ETz in the second 
frequency offset estimation processing becomes the range 
between -Rs/2Np and R^/lN^. 

Thus, it is understood that the frequency offset 
estimating range fnEn in the first frequency offset estimation 
processing is Np times larger than that in the second frequency 
offset estimation processing. That is, in the first frequency 
offset estimation processing, the frequency offset can be 
estimated in the relatively wide frequency offset estimating 
range fcETi- 

From the above description, in this Embodiment 1, both 
of the frequency f^ of the direct wave and the center frequency 
fji of the Doppler spread are estimated, and both of these 
frequencies f^ and f^ are used for compensation for the frequency 
offset, whereby deterioration in the BER characteristics can 
be prevented while securing a sufficient frequency offset 
compensating range. 

Hereinafter, the frequency offset estimation processing 
is explained in detail. The DSP 2 6 executes distortion amount 
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detection processing first (step Tl) . In the distortion amount 
detection processing, the distortion amounts of the 
transmission channel are detected in symbol units based on the 
known signals. That is, as shown by Formula (7) below, based 
on the Np-symbol known signal blocks included in the digital 
baseband signal r(kNF+i), the DSP 26 detects the distortion 
amounts GEpi(kNj,+i) of the transmission channel for each symbol 
in the known signal blocks. In this case, the distortion amounts 
Cj,pi(kNp+i) of the transmission channel corresponding to the 
amplitude and the distortion amount of the phase of the digital 
baseband signal r(kNp+i). In addition, in Formula (7) below, 
bp shows the symbol value of the known signal. 



Thereafter, the DSP 26 executes the first frequency 
offset estimation processing based on the detected distortion 
amounts CEpi(kNj,+i) of the transmission channel (step T2 ) . 
Concretely, the DSP 2 6 estimates the frequency of the direct 
wave as the first frequency offset based on the distortion 
amounts Cj,pi( kN^+i) of the transmission channel. More concretely, 
the DSP 26 estimates the phase difference eEpi(mLNF) 
corresponding to the phase rotation amount between the adjacent 
known signals as the first frequency offset based on the 
distortion amount CEpi(kNp+i) between the adjacent known signals 



CEpi(kNj. + i) = 



r(kNp + i) 
bp 



= Aa(kNp + i) + 



rXkNg. + i) 
bp 



(7) 
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among distortion amounts Cj,pi(kNp+i) of the transmission 
channel . 

Thereafter, the DSP 26 executes frequency offset 
elimination processing for eliminating the first frequency 
offset from the distortion amount of each known signal based 
on the estimated phase difference 9j,pi(mLNj.) (step T3 ) . More 
concretely, the DSP 2 6 compensates the phase by rotating the 
phases of the distortion amounts CEpi{kNp+i) of the transmission 
channel by the phase difference 0j,pi(mLNp) as shown by Formula 
(8). Thereby, the DSP 26 obtains the distortion amounts 
CEPi(kNj,+i) in which the first frequency offset has been 
eliminated. 

CEpi(kNp + i) = CEpi(kNp + i)exp[- jQ{kl\p + i)] 
e(kNp + i) = G(kNp + i - 1) + egj.i(mLNp) 

Thereafter, the DSP 2 6 executes the second frequency 
offset estimation processing (step T4 ) - Concretely, the DSP 26 
estimates the second frequency offset based on the distortion 
amounts CEpi(kN5,+i) in which the first frequency offset has been 
eliminated. More concretely, the DSP 26 estimates the phase 
difference 0EP2(mLNF) corresponding to the phase rotation amount 
for one symbol as the second frequency offset based on the 
distortion amount of the periodically inserted known signals 
among the distortion amounts CEpi(kNp+i) in which the first 
frequency offset has been eliminated. 

Thereafter, the DSP 2 6 executes phase synthesis 
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processing (step T5 ) . Concretely, the DSP 26 synthesizes the 
two phase differences eEpiCirLLNp) and 6Ep2(mLNp) as shown by Formula 
(9) to estimate the phase difference Gs(mLNp) as the last 
frequency offset. 

es(mLNF) = e^pi(mLNj.) + Q^^J,mL]<i^) (9) 

Fig. 7 is a flowchart for detailed explanation of the 
first frequency offset estimation processing. The DSP 26 
executes phase difference vector operation processing first 
(step Ul). Concretely, the DSP 2 6 determines the phase 
difference vector Dj.p(kNp) based on the distortion amounts 
Cj,pi(kNp+i) of the transmission channel. More concretely, as 
shown by Formula (10) below, based on the distortion amount 
between the adjacent known signals in an arbitrary known signal 
block among the distortion amounts Cj.pi(kNF+i) of the 
transmission channel, the DSP 26 determines the phase 
difference vector CEpi(kNp+i). In Formula (10) below, i shows 
the identification numbers corresponding one to one known 
signals, and are zero or larger and (Np-2) or smaller (0:5i^Np-2) . 
In addition, "*" shows a complex conjugate. 

Np-2 

DEp(kNj.) = 2^EP(i + l/^^F + ^ ^)^EPi * W + i) (10) 
i = 0 

Thereafter, the DSP 26 executes phase difference vector 
average processing (step U2 ) . Concretely, the DSP 26 averages 
the phase difference vectors DEp(kNp) in arbitrary known signal 
blocks for the estimation period T to determine the average 
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phase difference vector DEPA(inLNp) . 

More concretely, the DSP 26 collects the phase difference 
vectors DEp(kNp) in known signal blocks for the estimation period 
T from the estimation time (m-l)T to the estimation time mT (see 
Fig. 4) . If the number of known signal blocks in the estimation 
period T is L, the DSP 26 obtains L of phase difference vectors 
Dj,p{ (m-1 )T+lNp) . Herein, 1 shows the identification numbers 
corresponding one to one known signal blocks, and is 0 or more 
and (L-1) or less (O^lrsL-l). Thereafter, as shown by Formula 
(11) below, the DSP 26 determines the average phase difference 
vector DEpA(mLNp) by averaging the collected L of phase 
difference vectors Dep ( ( m- 1 ) LNp+lNp ) . 

L-1 

DEp^(mLNp) = ^D^^dm - l)LNp + IN^) (H) 

1=0 

Then, the DSP 2 6 executes phase difference operation 
processing (stepUS) . Concretely, as shown by Formula (12) below, 
the DSP 26 determines the phase difference 0Epi(mLNp) based on 
the average phase difference vector DEpA(mLNp) . Thus, the first 
frequency offset according to the frequency of the direct wave 
is estimated within the frequency offset estimating range fj^^^.^ 
that is wider than the frequency offset estimating range fDEr2' 

eWmLN,)=tan-^^^i5m(^^ (,2) 

Fig. 8 is a flowchart for explanation of the second 
frequency offset estimation processing. The DSP 26 executes 
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distortion amount average processing first (step VI). 
Concretely, the DSP 26 applies average processing to the 
distortion amounts Cj.pi(kNp+i) (O^i^Np-1) of the transmission 
channel in the known signal blocks in which the first frequency 
offset has been eliminated to determine the average distortion 
amount Cj.p ( kNp ) . 

More concretely, as shown by Formula (13) below, the DSP 
26 sums the distortion amounts CEpi(kNF+i) corresponding to each 
symbol in arbitrary known signal blocks, and divides the summed 
value by the symbol number Np of the known signals in the known 
signal blocks. Thus, the DSP 2 6 obtains the average distortion 
amount CEp(kNj,) of the transmission channel with regard to one 
known signal block in which noise or the like has been 
eliminated. 



CEp(kNp) = 




(13) 



Next, the DSP 26 executes interblock discrete Fourier 
transform (DFT) processing based on this average distortion 
amount CEp(kNj.) (step V2 ) . Concretely, the DSP 26 applies DFT 
processing to the average distortion amounts 
Cep( (m+l)Nj.) (O^l^L-l)in certain L of known signal blocks within 
the estimation period T (=LNj.Ts) , and thereby determines signal 
electric powers Pf(n) corresponding to the plurality of 
frequency offset candidates nAff^s. 

More concretely, the DSP 26 determines signal powers Pf (n) 
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corresponding to the plurality of frequency offset candidates 
nAfREs that are set for each predetermined estimation accuracy 
Afi^s as shown by the circles in Fig. 9A within the frequency 
offset estimating range fDEia* Herein, the frequency offset 
estimating range t^siir as mentioned above, is regulated by the 
known signal insertion period Nj,, and if the estimation 
accuracies Af^^ are used, this range is a range between -MAfRj,s 
and MAfRBs* M is a constant, and is approximately expressed by 
Formula (14) below. 

2NpAf^s ^ ^ 

The DSP 26 rotates the phases of the average distortion 
amounts c^p ( (m+1 )Nj,) in the frequency offset estimating range 
fDET2 by the phase amounts corresponding to the respective 
possible frequency offset candidates nAf^Es- Thereafter, the DSP 
2 6 vector-synthesizes the average distortion amounts 
Cep ( ( m+l ) Np ) whose phases have been rotated. Thereby, the DSP 
2 6 obtains signal powers P£(n) corresponding to the plurality 
of frequency offset candidates nAfj^g as shown by the group of 
arrows in Fig. 9B. 

The abovementioned processing can be summarized as 
Formula (15) below. In Formula (15) below, Rg shows the signal 
transmission rate, n is a value between -(M+W) and (M+W) . In 
addition, W is a parameter showing the frequency bandwidth of 
the frequency window described later. 
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2jilN^nAfj^ 



(15) 



Thereafter, the DSP 26 executes window signal power 
operation processing (step V3). Concretely, the DSP 26 
determines the window signal powers Ef(n) that correspond one 
to one frequency windows having predetermined frequency widths . 
For example, the frequency width of a frequency window is 2W 
times of the estimation accuracy Afj^s* W is set in accordance 
with the fading condition of the transmission channel, and is 
set to be, for example, two times of the Doppler spread. As shown 
by Formula (16) below, the DSP 26 sums the signal powers P£(n) 
of the frequency offset candidates nAfp^s existing in this 
frequency window to determine the window signal power Ef(n) 
corresponding to this frequency window. In Formula (16) below, 
n is a value between -M and M. 



Thus, by smoothing the signal powers Pf(n) by using the 
frequency windows, among the window signal powers Ef(n) 
corresponding to the plurality of frequency offset candidates 
nAfRj,g, as shown in Fig. 9C, the window signal power Ef(n) of 
the frequency offset candidate nAfRgg according to the center 
frequency of the Doppler spread becomes maximum. 

Next, the DSP 26 executes maximum value detection 
processing (stepV4) . Concretely, as shown by Formula (17) below. 




(16) 



k — w 
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the DSP 26 detects the maximum value among all window signal 
powers Ef (n) , and estimates the frequency offset candidate nAfj^g 
corresponding to the maximum value as the second frequency 
offset. That is, the DSP 26 estimates the center frequency of 
the Doppler spread as the second frequency offset. 



Next, the DSP 2 6 executes phase difference operation 
processing (stepVS). Concretely, as shown by Formula (18) below, 
the DSP 2 6 determines the phase difference 6^,52 (mLN^) of one 
symbol based on the estimated second frequency offset. Thus, 
the phase difference 9gp2(niLNp) corresponding to the center 
frequency of the Doppler spread is obtained. 



As described above, according to this Embodiment 1, by 
estimating both of the frequency of the direct wave and the 
center frequency of the Doppler spread, the last frequency 
offset is estimated. Therefore, a sufficient frequency offset 
compensating range can be secured, and excellent BER 
characteristics can be realized. That is, securing of a 
sufficient frequency offset compensating range and realization 
of excellent BER characteristics are compatible with each other. 
Therefore, even under a Rician fading environment, compensation 
for fading distortion can be satisfactorily carried out. 



Ef(nMAx) = M a x[Ef(n)] 



(17) 



9i 




(18) 
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Therefore, the original data can be restored with high quality. 
Accordingly, service for receiver users can be improved. 

Fig. 10 is a graph showing the frequency offset estimated 
characteristics relating to this Embodiment 1. In this figure, 
the circles show the BER in the case where the value obtained 
by synthesizing the frequency of the direct wave and the center 
frequency of the Doppler spread is regarded as frequency offset 
as in Embodiment 1. The triangles show the BER when only the 
frequency of the direct wave is regarded as frequency offset. 
Furthermore, the squares show the BER when only the center 
frequency of the Doppler spread is regarded as frequency offset. 

As clearly understood from this graph, the BER 
characteristics relating to Embodiment 1 are more excellent 
than those in the case where only the frequency of the direct 
wave is regarded as frequency offset, and a frequency offset 
compensating range, which is wider than that in the case where 
only the center frequency of the Doppler spread is regarded as 
frequency offset, is realized in Embodiment 1. 

Furthermore, according to the abovementioned Embodiment 
1 , frequency offset candidates nAfj^s are set for each estimation 
accuracy Af^j-s^ signal powers Pf(n) corresponding to the 
respective frequency offset candidates nAfj^s are determined, 
and the signal powers Pf (n) are smoothed by using the frequency 
windows, whereby the second frequency offset is estimated. 
Therefore, for example, if the estimation accuracies Af^Es are 
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finely divided, the estimation accuracy for the second 
frequency offset can be improved. Accordingly, further 
excellent BER characteristics can be realized. 

Embodiment 2 

Fig. 11 is a flowchart for explanation of the first 
frequency offset estimation processing relating to Embodiment 
2 of the invention. 

In the abovementioned Embodiment 1 , the phase difference 
eEpi(mLNp) is determined as the first frequency offset from the 
average phase difference vector DEp^(mLNy) . On the other hand, 
in this Embodiment 2, the phase difference 0^,51 (mLNp) is 
determined as the first frequency offset by further averaging 
the average phase difference vectors DEpa(mLNp) between the 
receiving burst signals by using a forgetting factor X, Thereby, 
the estimation accuracy for the frequency offset is improved. 

Concretely, the DSP 2 6 determines the phase difference 
vectors DEp(mLNF) (step Wl ) , averages the phase difference 
vectors DEp(mLNF) to determine the average phase difference 
vector DEp^CmLNj,) (step W2 ) , and then executes interburst average 
processing (step W3 ) . In the interburst average processing, the 
average phase difference vectors Dj,pA(niLNE,) are further averaged 
between the receiving burst signals. 

More concretely, the DSP 26 holds the average phase 
difference vector DEpa(mLNj,) at least until the next receiving 
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burst signal is received. Based on the average phase difference 
vector Dj,pA(mLNp) when the receiving burst signal B is received 
and the average phase difference vector <Depa ( mLNp ) >b_i when the 
one previous receiving burst signal (B-1) is received, the DSP 
26 determines the average phase difference vector <Depa { inLNp ) >b 
by way of Formula (19) below. In Formula (19) , X is a forgetting 
factor, and is a value of 0 or more and 1 or less (OsXsl). 
(D,,^(mLNp))^ = D,^JmLNp) + x(D^^JmLNp))^_^ (19) 
Thereafter, the DSP 26 executes phase difference 
operation processing based on the average phase difference 
vector <DEPA(mLNp)>B (step W4 ) . Concretely, as shown by Formula 
(20) below, the DSP 26 determines the phase difference 9Epi(mLNF) 
based on the average phase difference vector <Dj,pA ( inLNp ) >b . Thus, 
the DSP 26 estimates the first frequency offset coressponding 
to the frequency of the direct wave. 

, Im[(DEPA(niLI^))] 
Re[(DEpA(mLJ^))] 

As described above, according to this Embodiment 2, 
average processing is executed while the previous average phase 
difference vectors DEpj,(mLNp) are gradually forgotten by using 
the forgetting factor k. Therefore, even in the case where the 
fading condition changes with the times, average processing can 
be executed for the average phase difference vectors DEPA(nLLNp) 
while following the time variation. Accordingly, a phase 
difference vector in conformity with the condition of the 
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transmission path can be obtained. Therefore, even under an 
environment where the C/N is low and the frequency offset 
changes with the times, the frequency offset can be estimated 
with high accuracy. 

Embodiment 3 

Fig. 12 is a flowchart for explanation of the second 
frequency offset estimation processing relating to Embodiment 
3 of the invention. 

In the abovementioned Embodiment 1, the maximum value 
E£(nju^) of the signal power is detected among the window signal 
powers Ef(n) determined by one window signal power operation 
processing, whereby the second frequency offset is estimated. 
On the other hand, in this Embodiment 3, after averaging the 
window signal powers Ef(n) between the receiving burst signals 
by using the forgetting factor X, the maximum value is detected 
to estimate the second frequency offset. Thereby, the 
estimation accuracy of the frequency offset is improved. 

Concretely, the DSP 26 determines the average distortion 
amounts CEp(kNj,) (step XI), and applies interblock DFT 
processing to the average distortion amounts CEp(kNp), and 
thereby determines the signal powers Pf(n) (step X2 ) , and 
furthermore, determines the window signal powers Ef (n) from the 
signal powers Pf(n)(step X3 ) , and then executes window signal 
power average processing (step X4 ) . 
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More concretely, the DSP 2 6 holds the window signal power 
Ef (n) at least until the next receiving burst signal is received. 
Based on the window signal power Ef(n) when the receiving burst 
signal B is received and the average window signal power 
<Ef{n)>B_i when the one previous receiving burst signal (B-l) 
is received, the DSP 26 determines the average window signal 
power <Ef{n)>B by way of Formula (21) below. In Formula (21) 
below, X. is a forgetting factor between 0 and 1 (O^J^^l). 

<Ef(n))3 = E,(n) + K{-E,{n))^_^ (21) 

As shown by Formula (22) below, the DSP 26 detects the 
maximum value <Ef(n^x)>B average window signal power among the 
determined average window signal powers <E£(n)>B (step X5 ) . 

(Ef(nMAx))B = Max[(Ef(n))j (22) 

Thereafter, the DSP 26 estimates the phase difference 
6Ep2(mLNp) corresponding to the possible frequency offset nAfj^s 
corresponding to this maximum value <Ef(ntu^x)>B the second 
frequency offset (step X6). 

As described above, according to Embodiment 3, average 
processing is executed while the previous window signal powers 
are gradually forgotten by using the forgetting factor X,. 
Therefore, even in the case where the fading condition changes 
with the times, average processing of the window signal powers 
can be executed while following the time variation . Accordingly , 
window signal powers in conformity with the conditions of the 
transmission channel can be obtained. Therefore, even under an 
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environment where C/N is low and the frequency offset changes 
with the times, frequency offset can be estimated with high 
accuracy. 

Embodiment 4 

Fig. 13 is a flowchart for explanation of the second 
frequency offset estimation processing relating to Embodiment 
4 of the invention. 

In Embodiment 3, influences of noise are eliminated by 
averaging the window signal powers Ef (n) , whereby the estimation 
accuracy of the second frequency offset is improved. On the 
other hand, improvement in the estimation accuracy of the second 
frequency offset by eliminating influences of noise can also 
be achieved by averaging the signal powers P£(n). Therefore, 
in this Embodiment 4, the signal powers Pf(n) are averaged 
between the receiving burst signals to eliminate influences of 
noise, whereby the estimation accuracy of the second frequency 
offset is improved. 

Concretely, the DSP 26 determines the average distortion 
amounts CEPACkNp) (step Yl) and applies DFT processing to the 
average distortion amounts Cj,pA(kNp) to determine signal powers 
PjE(n) (step Y2 ) , and then executes signal power average 
processing (step Y3 ) . More concretely, the DSP 26 holds the 
signal power Pf (n) at least until the next receiving burst signal 
is received. Based on the signal power P£(n) when the receiving 
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burst signal B is received and the average signal power <Pf {n)>B_i 
when the one previous receiving burst signal (B-1) is received, 
the DSP 26 determines the average signal power <Pf(n)>B by way 
of Formula (23) below. In Formula (23) , X Is a forgetting factor 
between 0 and 1 (O^eXsl). 



Thereafter, the DSP 26 determines the window signal 
powers Ef(n) based on average signal powers <Pf(n)>B thus 
determined ( step Y4 ) , and then detects the maximum value ( n^^^) 
window signal power (step Y5), and estimates the phase 
difference 6EP2(mLNp) corresponding to the frequency offset 
candidate nAfj^s corresponding to the maximum value Ef(n^^^) as 
the second frequency offset (step Y6). 

As described above, according to this Embodiment 4, the 
average processing is executed while the previous signal powers 
Pf(n) are gradually forgotten by using the forgetting factor 
X. Therefore, even in the case where the fading condition 
changes with the times, the average processing of the signal 
powers can be executed while following the time variation. 
Accordingly signal powers in conformity with the conditions of 
the transmission channel can be obtained. Therefore, even under 
an environment where C/N is low and the frequency offset changes 
with the times, the frequency offset can be estimated with high 
accuracy. 




(23) 
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Embodiment 5 

Fig. 14 is a flowchart for explanation of the first 
frequency offset estimation processing relating to Embodiment 
5 of the invention. 

In the abovementioned Embodiment 1, based on the 
distortion amounts CEpi(kNp+i) of the transmission channel, the 
phase difference vector DepCIcNf) between the adjacent known 
signals is determined. On the other hand, in this Embodiment 
5 by deteraining only phase difference information that is the 
scalar quantity between the adjacent known signals based on the 
distortion amounts Cj.pi ( kNj,+i ) of the transmission channel, the 
processing is simplified. 

Concretely, the DSP 2 6 executes phase difference 
information operation processing first between the known 
signals (step Zl). In this phase difference information 
operation processing, phase difference information as the 
scalar quantity between adjacent known signals is determined. 
That is, as shown by Formula (24) below, the DSP 26 converts 
the distortion amounts Cjpi(kNp+i) of the transmission channel 
into phase information 6EPi(kNp+i), and determines phase 
difference information AeEp(kNj,) between adjacent known signals . 
In Formula (24) below, i shows the identification numbers which 
correspond one to one known signals, and is a value between 0 
and (Np-2) (O^i^Np-2). 
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1 

AeEp(kN^) = 2(eEPi(kNp + i + 1) - eEpi(kNp + i)) 

'^p ~ ■•■ 1=0 

e,,,(kN, + i) = tan- ^"^[^KPi(kN, .i)] (24) 
Re[CEpi(kNj. + i) ] 

Thereafter, the DSP 26 executes average processing for 

averaging the phase difference information AeBp( kNj,) in one burst 

(step Z2). Concretely, the DSP 26 collects phase difference 

information Aej.p(kNp) for the estimation period T from the 

estimation time (m-l)T to the estimation time mT. If the number 

of known signal blocks in the estimation period T is L, the DSP 

26 obtains L of phase difference information Aej,p( ( m- 1 ) LN^+lNp ) . 

Herein, 1 shows the indication numbers which correspond one to 

one known signal blocks, and is a value between 0 and (L-1) 

(O^lsL-1) . 

Thereafter, as shown by Formula (25) below, the DSP 26 
applies average processing to the collected L of phase 
difference information Ae^p ( (m-1 )LNj,+lNj,) , whereby the average 
phase difference egpi(mLNp) is determined as the first frequency 
offset. 

2 L-l 

eEPi(mLNp) = - 2) ^9ep( (m - l)LNp + IN^.) (25) 

^ 1=0 

As described above, according to this Embodiment 5, the 
average processing is executed by using only the phase 
difference information. Therefore, the processing can be 
simplified. 
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Embodiment 6 

Fig. 15 is a flowchart for explanation of the first 
frequency offset estimation processing relating to Embodiment 

6 of the invention. 

In the abovementioned Embodiment 1, the distortion 
amounts CEpi(kNp+i) of the transmission channel are multiplied 
by the complex conjugate and summed, whereby the phase 
difference vector D^pCkNp) between adjacent known signals is 
determined. On the other hand, in Embodiment 6, by applying DFT 
processing to the distortion amounts Cj,pi( kNj.+i) of the 
transmission channel, the phase difference vector Dj,p(kNp) is 
determined . 

Concretely, the DSP 2 6 executes inter-symbol DFT 
processing (step Rl ) . More concretely, by applying the DFT 
processing to the distortion amounts CEpi(kNp+i) of the 
transmission channel, the DSP 26 determines signal electric 
powers Pfi(n) corresponding to the plurality of frequency offset 
candidates nAfj^s. 

Still more concretely, as shown by Formula (26), within 
the frequency offset estimating range f^.^^, between -Rs/2 and Rs/2 , 
the DSP 26 rotates the phases of the distortion amounts 
Cj,pi( (m-1 )LNp+lNj, +i) of the transmission channel by the phase 
amounts corresponding to the respective frequency offset 
candidates nAf^^g. Thereafter, the DSP 26 vector-synthesizes the 
distortion amounts CEpi( (m-1 )LNp+lNp +i) whose phases have been 
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rotated, and determines signal powers Pfi(n) corresponding to 
the plurality of frequency offset candidates nAf^Es within the 
frequency of f set estimating range. In addition, in Formula (26) 
below, n is approximately a value between -R^/ {2Af^s) and 
R,/(2Af^3). 

P,,(n) = 2I 2'c,p,((m - 1)LN, + IN, + i)expf- j ^"^^^^ ^ ^^^^^ 

i=o|i=o y Rg 

(26) 

Then, the DSP 26 vector-synthesizes the distortion 
amounts c^^^{ (m-1 )LNj.+lNj. +i) whose phases have been rotated, and 
obtains signal powers Pfi(n) corresponding to the respective 
plurality of frequency offset candidates nAfj^s- 

Then, the DSP 26 executes window signal power operation 
processing (step R2 ) . Concretely, as shown by Formula (27) , the 
DSP 2 6 sums the signal powers Pfi(n) corresponding to the 
frequency offset candidates nAf^Es in the window having a 
frequency width of 2WiAfj^s, whereby the window signal powers 
Efi(n) are determined. 

Efi(n) = + (27) 

Thereafter, the DSP 26 executes maximum value detection 
processing to detect the maximum value E^^{n^) among the window 
signal powers E^^{n) (step R3 ) , and executes phase difference 
operation processing based on the detected maximum value 
Efi(nj^) (step R4), whereby the phase difference 0j.pi(mLNF) is 
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estimated as the first frequency offset. 

As described above, according to Embodiment 6, frequency 
offset candidates nAfuj-g are set for each estimation accuracy 
AfRj,s, signal powers Pfi(n) corresponding to the respective 
frequency offset candidates nAfj^g are determined, and the 
signal powers Pfi(n) are smoothed by using the frequency window, 
whereby the first frequency offset is estimated. Therefore, for 
example, if the estimation accuracies AfR^s are finely set, the 
estimation accuracy of the first frequency offset can be 
improved . 

Embodiment 7 

Fig. 16 is a flowchart for explanation of the frequency 
offset estimation processing relating to Embodiment 7 of the 
invention. 

In the abovementioned Embodiments 1 through 6, the 
distortion amounts Cj,pi(kNj.+i) of the transmission channel are 
compensated based on the phase difference 9j,pi(mLNp) as the first 
frequency offset, and the phase difference 9EP2(niLNj,) as the 
second frequency offset based on the distortion amounts 
Cj.pi(kNp+i) after being compensated is estimated. 

Herein, the first frequency offset corresponds to the 
frequency f^ of the direct wave, and the frequency f^ of the 
direct wave exists near the center frequency f„ of the Doppler 
spread (see Fig. 9A) . Therefore, on supposition that a certain 
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range around the frequency of the direct wave is a frequency 
offset estimating range, it can be expected that the center 
frequency f^ of the Doppler spread exists within the range. 

Therefore, in Embodiment 7, the frequency offset 
estimating range for estimation of the second frequency offset 
based on the first frequency offset is regulated. 

Concretely, the DSP 26 detects the distortion amounts 
CEpi(kNj,+i) of the transmission channel (step Ql), and then 
executes the first frequency offset estimation processing and 
the second frequency offset estimation processing 
simultaneously. More concretely, the DSP 26 estimates the first 
frequency offset based on the abovementioned distortion amounts 
CEpi(kNj,+i) of the transmission channel. The DSP 26 temporarily 
holds, for example, this first frequency offset. 

In addition, the DSP 26 executes the first frequency 
offset estimation processing and distortion amount average 
processing simultaneously (step Q2 ) . Concretely, the DSP 26 
applies average processing to the distortion amounts CEpi(k;Np+i) 
of the transmission channel to determine the average distortion 
amount Cj.p ( kNp ) . 

Thereafter, the DSP 26 executes DFT processing (steps Q3 
and Q4 ) . Concretely, the DSP 26 applies DFT processing to L of 
average distortion amounts c^^{ (m+l)Np) in the estimation period 
T (=LNj.Ts), and determines signal powers Pf(n) corresponding to 
the plurality of frequency offset candidates nAfj^s 
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respectively . 

More concretely, the DSP 26 sets a frequency offset 
estimating range based on the first frequency offset (step Q3 ) , 
and determines signal powers Pf(n) corresponding to the 
plurality of possible frequency offsets nAfj^g within the 
frequency offset estimating range respectively (step Q4 ) . 

Still more concretely, the DSP 26 accords the first 
frequency offset to any one of the frequency offset candidates 
nAfREs set in advance for each estimation accuracy AfR^s* Herein, 
the frequency offset candidate in accordance with the first 
frequency offset is shown by NAfjjEg. In addition, the DSP 26 sets 
the range around this frequency offset candidate NAfj^s/ that 
is, the range between (-M+N)AfRj,s and (M+N)AfRj.s as the frequency 
offset estimation range. M is the constant shown in Formula 



Thereafter, as shown by Formula (28) below, based on the 
phase amounts corresponding to the respective frequency offset 
candidates Afj^g within this set frequency offset estimation 
range, the DSP 26 rotates the phases of the average distortion 
amounts Cj.p ( ( m+1 ) Nj. ) (OslsL-l). In Formula (28) below, n is a 
value between (-M-W+N) and (M+W+N) . 



The DSP 26 vector-synthesizes the average distortion 
amounts Cj,p( (m+1 )Nj,) whose phases have been rotated. Thereby, 



(14). 
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the DSP 2 6 determines signal powers Pf(n) corresponding to the 
plurality of frequency offset candidates nAfj^g respectively. 

Thereafter, the DSP 26 determines window signal power 
Ef(n) by cyclically summing the signal powers P£(n) in a 
frequency window with a predetermined frequency width ( step Q5 ) . 
The DSP 26 determines the maximum value Ef(nHjix) of the window 
signal power (step Q6), and estimates the frequency offset 
candidate nAfj^g of this maximum value EfCnj^x) as the last 
frequency offset (step Q7 ) . 

As described above, according to Embodiment 1, the 
frequency offset estimating range to estimate the frequency 
offset based on the first frequency offset is set. Therefore, 
the first frequency offset elimination processing and phase 
synthesis processing become unnecessary. Therefore, the 
processing can be simplified. 

In addition, since the first frequency offset is 
estimated in a relatively wide frequency offset estimating 
range, a sufficient frequency offset compensating range is 
secured. Also, since the center frequency of the Doppler spread 
is used to estimate the last frequency offset, excellent BER 
characteristics are secured. Therefore, the same effects as in 
Embodiment 1 such as securing of a sufficient frequency offset 
compensating range and prevention of deterioration in the BER 
characteristics can also be obtained. 
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Embodiment 8 

Fig. 17 is a flowchart for explanation of demodulation 
processing relating to Embodiment 8 of the invention. 

In Embodiment 1, the frequency offset elimination 
processing in which the phase of the digital baseband signal 
is rotated is executed after filtering processing. On the other 
hand, in this Embodiment 8, the frequency offset elimination 
processing is executed previous to the filtering processing. 

Frequency offset of the receiving burst signal occurs due 
to the stability of the oscillation circuit of the transmitter 
1 and receiver 10 as explained in the subsection "Discription 
of the Related Art" above. Particularly, in the case where 
low-cost oscillation circuits are used, since the frequency 
stability is low, great frequency offset occurs. In addition, 
the frequency stability of the oscillation circuits is 
influenced by the surrounding environment such as temperature 
and changes in the power supply voltage. 

Therefore, there is every possibility that the frequency 
offset becomes greater in accordance with time elapse although 
the offset is slight at the beginning. In this case, the 
frequency bandwidth of the receiving burst signal may expand 
to a large frequency bandwidth in comparison with the cutoff 
frequency in the filtering processing. Therefore, a part of the 
receiving burst signal may be cut by the filtering processing. 
Accordingly, the data decision may not be sufficiently carried 
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out. 

In view of the above circumstances, in this Embodiment 
8, the frequency offset elimination processing is executed 
prior to the filtering processing. Concretely, in this 
Embodiment 8, the frequency offset elimination processing is 
executed prior to the filtering processing, and distortion 
amount detection processing, frequency offset estimation 
processing, and integration processing are executed 
simultaneously with sequential processing including fading 
distortion compensation processing and data decision 
processing. 

More concretely, as shown in Fig. 17A, when digital 
baseband signals are supplied from the A/D converter circuits 
25a and 25b, previous to the filtering processing, the DSP 26 
executes frequency offset elimination processing (step Pi). 
Still more concretely, based on the accumulated phase 
difference e( (m-1 )LNp+lNF+i) determined by the separately 
executed integration processing, the phases of the digital 
baseband signals are rotated in reverse. Thereby, digital 
baseband signals from which frequency offset has been 
eliminated are obtained. 

Thereafter, the DSP 26 applies filtering processing to 
the digital baseband signals (step P2 ) to eliminate noise 
components and the like. Thereafter, the DSP 2 6 detects a 
digital baseband signal corresponding to the Nyquist point 
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(step S3), eliminates fading distortion from the digital 
baseband signal {step P4 ) , and then executes data decision 
processing (step P5). 

On the other hand, as shown in Fig. 17B, the DSP 26 decides 
whether or not the digital baseband signal corresponding to the 
Nyquist point has been detected (step Nl ) . If the Nyquist point 
has been detected, based on the digital baseband signal 
corresponding to the Nyquist point, the DSP 26 detects the 
distortion amount Cj.pi(kNp+i) of the transmission channel (step 
N2). Thereafter, based on the distortion amount Cj.pi(kNp+i) of 
the transmission channel, the DSP 26 estimates the phase 
difference 9s(mLNj,) corresponding to the frequency offset (step 
N3 ) , and cyclically sums the phase differences es(mLNp,) , whereby 
an accumulated phase dif ference e( (m-1 )LNp+lNp+i) is determined 
(step N4). The DSP 26 uses this determined accumulated phase 
difference 6 ( (m-1 )LNp+lNF+i) for the frequency offset 
elimination processing in step Pi. 

As described above, according to this Embodiment 8, the 
frequency offset elimination processing is executed before the 
filtering processing. Therefore, even if the frequency 
bandwidths of the receiving burst signals are greater than the 
cutoff frequency in the filtering processing, the frequency 
offset can be eliminated while a part of the receiving burst 
signals is not cut. Therefore, data decision can be 
satisfactorily carried out. 
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Embodiment 9 

Fig. 18 is a flowchart for explanation of demodulation 
processing relating to Embodiment 9 of the invention. 

In the abovementioned Embodiments 1 through 8, by 
directly rotating the phase of the digital baseband signal, the 
frequency offset is compensated. On the other hand, in this 
Embodiment 9, the phase of the receiving burst signal is rotated 
by changing the frequency of the local oscillation signal 
generated at the oscillation circuit 22, whereby the frequency 
offset is compensated. 

Concretely, in this Embodiment 9, digital signal 
processing by way of software is not applied to the digital 
baseband signal, but a voltage to be applied to the oscillation 
circuit 22 comprised of VCO is controlled to eliminate the 
frequency offset. 

That is, in this Embodiment 9, as shown in Fig. 18A, 
filtering processing (step Ml), Nyquist point detection 
processing (step M2), fading distortion compensation 
processing (step M3), and data decision processing (step M4 ) 
are executed, and simultaneously with these sequential 
processings, as shown in Fig. 18B, distortion amount detection 
processing, frequency offset estimation processing, and VCO 
control processing for controlling the voltage to be applied 
to the VCO 22 are executed. 
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More concretely, the DSP 26 decides whether or not the 
digital baseband signal corresponding to the Nyquist point has 
been detected (step LI). In the case where the Nyquist point 
has been detected, the DSP 26 detects the distortion amount 
CEpi(kNj,+i) of the transmission channel based on the digital 
baseband signal corresponding to the Nyquist point (step L2). 
Thereafter, the DSP 26 estimates frequency offset based on the 
distortion amount CEpi(kNp+i) of the transmission channel (step 
L3) . Then, the DSP 26 supplies a voltage in accordance with the 
estimated frequency offset to the VCO 22 (step L4 ) . 

Thereby, the VCO 22 oscillates a local oscillation signal 
having an oscillation frequency which has deviated by an amount 
according to the frequency offset from the original oscillation 
frequency. Therefore, the frequency converter circuit 21 
outputs an analog baseband signal from which frequency offset 
has been eliminated. Thus, excellent frequency offset 
compensation can be realized. 

As described above, according to Embodiment 9, by 
adjusting the frequency of the local oscillation signal 
generated by the VCO 22 , the frequency offset is eliminated from 
the receiving IF signal. That is, before the filtering 
processing by the DSP 26, the frequency offset is eliminated. 
Therefore, as in the abovementioned Embodiment 8, even if the 
frequency bandwidth of the receiving burst signal expands to 
be larger than the cutoff frequency in the filtering processing. 
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frequency offset can be eliminated while a part of the receiving 
burst signal is not cut. Therefore, the data decision can be 
satisfactorily carried out. 

Embodiment 10 

Fig. 19 is a flowchart for explanation of demodulation 
processing relating to Embodiment 10 of the invention. 

In the abovementioned Embodiment 1, the first frequency 
offset and the second frequency offset are synthesized to 
estimate the last frequency offset, and in accordance with this 
frequency offset, compensation for the frequency offset is 
carried out, whereby the frequency of the digital baseband 
signal is automatically controlled. On the other hand, in this 
Embodiment 10, frequency offset compensation in accordance with 
the first frequency offset is carried out for the digital 
baseband signals, and frequency offset compensation in 
accordance with the second frequency offset is carried out for 
digital baseband signals after the abovementioned compensation 
is carried out, whereby the frequency of the digital baseband 
signal is automatically controlled. 

Concretely, the DSP 2 6 applies predetermined filtering 
processing to the digital baseband signals (step Kl ) , and then 
detects a digital baseband signal corresponding to Nyquist 
point (step K2 ) . Next, the DSP 2 6 executes the first automatic 
frequency control processing (hereinafter, referred to as the 
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first AFC processing) for the digital baseband signal (stepKS), 
and then executes the second automatic frequency control 
processing (hereinafter, referred to as the second AFC 
processing) for the digital baseband signal after being 
subjected to the first AFC processing (step K4 ) . Then, the DSP 
26 executes fading distortion compensation processing for the 
digital baseband signal after being subjected to the second AFC 
processing (step K5), and applies data decision processing to 
the signal ( step K6 ) , whereby data corresponding to the original 
information signal is restored. 

Fig. 20 is a flowchart for explanation of the first AFC 
processing. The DSP 26 detects the distortion amounts Cj,i,i( kNp+i) 
of the transmission channel for each symbol (step Jl ) , and then 
executes the first frequency offset estimation processing (step 
J2). The first frequency offset estimation processing is the 
same as that in the abovementioned Embodiment 1, and by this 
processing, the DSP 26 estimates the phase difference eEpi(mLNp) 
as the first frequency offset corresponding to the phase 
rotation amount for one symbol. 

Thereafter, the DSP 26 executes integration processing 
for cyclically summing the phase differences eBpi(mLNp) (step 
J3 ) to obtain an accumulated phase difference (m-1 )LNF+lNj.+i) . 
In this case, the DSP 26 obtains the phase difference 
6i( (m-1 )LNF+lNp+i) accumulated from the estimation time mT for 
each symbol period Tg. 
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Thereafter, the DSP 26 executes the first frequency 
offset elimination processing (step J4 ) . Concretely, based on 
the accumulated phase difference (m-1 )LNy+lNj.+i) , the DSP 26 
rotates in reverse the phase of the digital baseband signal 
r{kNp+i) by an amount in accordance with the accumulated phase 
difference e,( (m-1 )LNp+lNp+i) . Thereby, the first frequency 
offset can be eliminated from the digital baseband signal 
r{kNp+i) . Thus, the DSP 2 6 achieves automatic frequency control 
of the digital baseband signals r(kNp+i) based on the first 
frequency offset. 

Fig. 21 is a flowchart for explanation of the second AFC 
processing. Based on the known signal blocks of the digital 
baseband signal after being subjected to the first AFC 
processing, the DSP 26 detects the distortion amount CEpi2(kNF+i) 
of the transmission channel for each symbol in the known signal 
blocks (step II). 

Thereafter, based on the detected distortion amount 
Cj,pi2(kN5,+i) of the transmission channel, the DSP 2 6 executes the 
second frequency offset estimation processing that is the same 
as in Embodiment 1 (step 12) . Thereby, the DSP 2 6 estimates the 
phase difference eEP2(mLNj.) as the second frequency offset 
corresponding to the phase rotation amount for one symbol. 

Thereafter, the DSP 26 executes integration processing 
for cyclically summing the phase differences ej,p2(mLNp) (step 
13), and obtains an accumulated phase difference e2((m- 
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1 )LNp+lNp+i) . In this case, the DSP 26 obtains the phase 
difference ( (m-1 )LNp+lNj.+i) accumulated from the estimation 
time mT for each symbol period Tg. 

Thereafter, the DSP 26 executes the second frequency 
offset elimination processing (step 14). Concretely, based on 
the accumulated phase difference Q^i (m-1 )LN5,+lN5,+i) , the DSP 26 
rotates in reverse the phase of the digital baseband signal 
rRi(kNp+i) by an amount according to the accumulated phase 
difference (m-1 )LNp+lNF+i) . Thereby, the second frequency 
offset can be eliminated from the digital baseband signal 
rRi(kNp+i). Thus, the DSP 26 achieves automatic frequency 
control of the digital baseband signal rRi(kNp+i) based on the 
second frequency offset. 

As described above, according to this Embodiment 10, 
compatibility of securing of a sufficient frequency offset 
compensating range and securing of excellent BER 
characteristics is also achieved. 

Embodiment 11 

Fig. 22 is a flowchart for explanation of demodulation 
processing relating to Embodiment 11 of the invention. 

In the abovementioned Embodiment 10 , the frequency offset 
elimination processing is executed after the filtering 
processing. On the other hand, in this Embodiment 11, the 
frequency offset elimination processing is executed before the 
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filtering processing. 

Concretely, in this Embodiment 11, the first frequency 
offset elimination processing that is a part of the first AFC 
processing is executed before the filtering processing, and the 
remaining processing of the first AFC processing, that is, 
distortion amount detection processing, the first frequency 
offset estimation processing, and integration processing are 
executed simultaneously with sequential processing. 

More concretely, as shown in Fig. 22A, when the digital 
baseband signals are supplied from the A/D converter circuits 
25a and 25b, and previous to the filtering processing, the DSP 
26 executes the first frequency offset elimination processing 
(step HI) . Still more concretely, the DSP 26 rotates in reverse 
the phases of the digital baseband signals by amounts according 
to the accumulated phase difference e( (m-1 )LNp+lNj.+i) 
corresponding to the first frequency offset determined by the 
integration processing executed simultaneously with the 
sequential processing. Thereby, digital baseband signals from 
which first frequency offset has been eliminated are obtained. 

Thereafter, the DSP 26 applies filtering processing to 
the digital baseband signals (step H2 ) to eliminate noise 
components and the like, and then detects a digital baseband 
signal corresponding to the Nyquist point ( step H3 ) . Thereafter, 
the DSP 26 executes the second AFC processing based on the 
digital baseband signal (step H4 ) . By executing this second AFC 
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processing, the DSP 26 can obtain a digital baseband signal in 
which the last frequency offset has been eliminated. Next, the 
DSP 26 eliminates fading distortion from the digital baseband 
signal (step H5), and then executes data decision processing 
(step H6) . 

On the other hand, as shown in Fig. 22B, the DSP 26 decides 
whether or not the digital baseband signal corresponding to the 
Nyquist point has been detected (step Gl ) . In the case where 
the Nyquist point has been detected, based on the digital 
baseband signal corresponding to the Nyquist point, the DSP 26 
detects the distortion amount CEpi(kNj,+i) of the transmission 
path (step G2). Thereafter, based on the distortion amount 
CEpi(kNj.+i) of the transmission channel, the DSP 26 estimates 
the phase difference es(mLN5,) corresponding to the first 
frequency offset (step G3 ) , and cyclically sums the phase 
differences es(mLNj,), whereby an accumulated phase difference 
e( (m-l)LNF+lNp+i) is determined (step G4 ) . The DSP 26 uses this 
determined accumulated phase difference 6 ( (m-1 )LNp+lNp+i) for 
the first frequency offset elimination processing in step Hi. 

As described above, according to this Embodiment 11, the 
first frequency offset elimination processing is executed 
before the filtering processing. Therefore, even when the 
frequency bandwidth of the receiving IF signal becomes larger 
than the cutoff frequency in the filtering processing, the first 
frequency offset can be eliminated while a part of the receiving 
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IF signals is not cut. Therefore, data decision can be 
satisfactorily carried out. 

Embodiment 12 

Fig. 23 is a flowchart for explanation of demodulation 
processing relating to Embodiment 12 of the invention. 

In the abovementioned Embodiment 10 , by directly rotating 
the phase of the digital baseband signal, the frequency offset 
is compensated. On the other hand, in this Embodiment 12, by 
adjusting the frequency of the local oscillation signal 
generated by the oscillation circuit 22 that is the VCO, the 
phase of the receiving burst signal is rotated, whereby the 
frequency offset is compensated. 

Concretely, in this Embodiment 12, digital signal 
processing by way of software is not applied to the digital 
baseband signal, but the voltage to be applied to the 
oscillation circuit 22 formed of the VCO is controlled, whereby 
the frequency offset is eliminated. 

That is, in this Embodiment 12, as shown in Fig. 23A, the 
filtering processing (step Fl), Nyquist point detection 
processing (step F2 ) , the second AFC processing (step F3 ) , 
fading distortion compensation processing (step F4 ) , and data 
decision processing (step F5) are executed, and simultaneously 
with these sequential processings, as shown in Fig. 23B, 
distortion amount detection processing, the first frequency 
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offset estimation processing, and VCO control processing for 
controlling the voltage to be applied to the VCO 22 are executed. 

Concretely, the DSP 26 decides whether or not the digital 
baseband signal corresponding to the Nyquist point has been 
detected (step El ) . In the case where the Nyquist point has been 
detected, based on the digital baseband signal corresponding 
to the Nyquist point, the DSP 26 detects the distortion amount 
CEpi{kNj,+i) of the transmission channel (step E2 ) . Thereafter, 
the DSP 26 estimates the first frequency offset based on the 
distortion amount Cj.pi ( kNp+i ) of the transmission channel (step 
E3). Then, the DSP 26 applies a voltage according to the 
estimated first frequency offset to the VCO 22 (step E4). 

Thereby, the VCO 22 oscillates the local oscillation 
signal having an oscillation frequency deviated by an amount 
in accordance with the first frequency offset from the original 
oscillation frequency. Therefore, the frequency converter 
circuit 21 outputs analog baseband signals in which the first 
frequency offset has been eliminated. Thus, excellent frequency 
offset compensation can be realized. 

As described above, according to this Embodiment 12, by 
adjusting the frequency of the local oscillation signal 
generated by the VCO 22, frequency offset is eliminated from 
the receiving IF signal. That is, the first frequency offset 
is eliminated before the filtering processing. Therefore, as 
in the abovementioned Embodiment 1 1 , even when the frequency 
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bandwidth of the receiving IF signal expands to be larger than 
the cutoff frequency in the filtering processing, frequency 
offset can be eliminated while a part of the receiving IF signal 
is not cut. Therefore, data decision can be satisfactorily 
carried out. 

Other Embodiments 

The Embodiments of the invention are as described above, 
however, the invention is not limited to the abovementioned 
Embodiments. For example, in the Embodiments, the case where 
the TDMA is applied as the communications method is described. 
However, as the communications method, the FDMA (Frequency 
Division Multiple Access) , CDMA (Code Division Multiple Access ) , 
or the like may be applied. In such a case, the radio receiving 
signals are not burst signals but are continuous signals as in 
the case where TDMA is applied. However, needless to say, the 
invention can be easily applied even to this case. 

In addition, in the abovementioned Embodiments, the case 
where the demodulation processing is realized by way of software 
by the DSP 2 6 is described. However, of course, the steps of 
the demodulation processing to be executed by the DSP 26 may 
be realized by hardware circuits. 
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WHAT IS CLAIMED IS: 

1. An automatic frequency control method for controlling 
the frequency of a radio receiving signal by compensating 
frequency offset of the radio receiving signal periodically 
including a plurality of adjacent known signals, wherein the 
frequency of the direct wave of the radio receiving signal and 
the center frequency of the Doppler spread of the radio 
receiving signal are estimated based on the distortion amounts 
of the known signals included in the radio receiving signal, 
and the frequency offset of the radio receiving signal is 
compensated based on both of these frequencies. 

2. An automatic frequency control method according to 
claim 1, wherein the radio receiving signal is a burst signal 
synchronized with a predetermined time slot in the TDMA. 

3 . An automatic frequency control method for controlling 
the frequency of a radio receiving signal by eliminating 
frequency offset from the radio receiving signal periodically 
including a plurality of adjacent known signals, comprising: 

a frequency offset estimating step for estimating 
frequency offset of the radio receiving signal from the 
frequency of the direct wave of the radio receiving signal, 
which is estimated based on the distortion amounts of the known 
signals included in the radio receiving signal, and the center 
frequency of the Doppler spread of the radio receiving signal; 
and 
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a frequency offset eliminating step for eliminating the 
estimated frequency offset from the radio receiving signal. 

4 . An automatic frequency control method according to 
claim 3, wherein the frequency offset estimating step includes: 

a distortion amount operating step for determining the 
distortion amounts of the known signals included in the radio 
receiving signal; 

a first frequency offset estimating step for estimating 
the frequency of the direct wave of the radio receiving signal 
based on the distortion amounts between adjacent known signals 
among these determined distortion amounts; 

a distortion amount frequency offset eliminating step for 
eliminating frequency offset from the distortion amounts of the 
known signals based on the estimated frequency of the direct 
wave; 

a second frequency offset estimating step for estimating 
the center frequency of the Doppler spread of the radio 
receiving signal based on the distortion amounts between the 
periodically inserted known signals among the distortion 
amounts from which the frequency offset has been eliminated; 
and 

a synthesizing step for synthesizing the phase amount 
corresponding to the frequency of the direct wave and the phase 
amount corresponding to the center frequency of the Doppler 
broadening to estimate frequency offset of the radio receiving 
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signal . 

5. An automatic frequency control method according to 
claim 4, wherein the first frequency offset estimating step 
includes steps of: 

determining the phase difference vectors based on the 
distortion amounts between adjacent known signals among the 
determined distortion amounts ; 

determining the average phase difference vector by 
averaging the determined phase difference vectors for a 
predetermined period; and 

estimating the frequency of the direct wave based on the 
determined average phase difference vector. 

6 . An automatic frequency control method according to 
claim 4, wherein the second frequency offset estimating step 
includes steps of: 

determining the average distortion amount by averaging 
the distortion amounts from which frequency offset has been 
eliminated; 

determining signal powers corresponding to a plurality 
of frequency offset candidates set at each predetermined 
interval within a predetermined frequency offset estimating 
range based on the determined average distortion amount; 

determining the window signal power of all frequency 
offset candidates within the frequency offset estimating range 
by summing the signal powers of the frequency offset candidates 
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in the frequency window with a predetermined frequency width 
among the determined signal powers; and 

estimating the frequency offset candidate corresponding 
to the maximum value of the determined window signal power as 
the center frequency of the Doppler spread. 

7 . An automatic frequency control method according to 
claim 4, wherein the first frequency offset estimating step 
includes steps of: 

determining phase difference information based on the 
distortion amounts between the adjacent known signals among the 
distortion amounts; 

determining average phase difference information by 
averaging the determined phase difference information for a 
predetermined period; and 

estimating the frequency of the direct wave based on the 
determined average phase difference information. 

8- An automatic frequency control method according to 
claim Af wherein the first frequency offset estimating step 
includes steps of: 

determining signal powers corresponding to a plurality 
of frequency offset candidates set at each predetermined 
interval in a predetermined frequency offset estimating range 
based on the distortion amounts of the adjacent known signals 
among the distortion amounts; 

determining window signal powers corresponding to all 
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frequency offset candidates by summing signal powers 
corresponding to the frequency offset candidates in a frequency 
windows with predetermined frequency widths among the 
determined signal powers; and 

estimating the frequency offset candidate corresponding 
to the maximum value of the determined window signal powers as 
the frequency of the direct wave. 

9. An automatic frequency control method according to 
claim 3, wherein the frequency offset estimating step includes: 

a distortion amount operating step for determining 
distortion amounts of the known signals included in the radio 
receiving signal; 

a first frequency offset estimating step for estimating 
the frequency of the direct wave of the radio receiving signal 
based on the distortion amounts between the adjacent known 
signals among the determined distortion amounts; 

an average distortion amount operating step for 
determining the average distortion amount by averaging the 
determined distortion amounts; 

a signal power operating step for determining signal 
powers corresponding to a plurality of frequency offset 
candidates set at each predetermined interval in a frequency 
offset estimating range regulated by the estimated frequency 
of the direct wave based on the determined average distortion 
amount ; 
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a window signal power operating step for determining 
window signal powers of all frequency offset candidates in the 
frequency offset estimating range by summing the signal powers 
of the frequency offset candidates in the frequency window with 
a predetermined frequency width among the determined signal 
powers ; and 

a second frequency offset estimating step for estimating 
the frequency offset candidate corresponding to the maximum 
value of the determined window singal powers as the frequency 
offset of the radio receiving signal. 

10. An automatic frequency control method according to 
claim 3, further comprising a filtering step for eliminating 
high frequency components over the cutoff frequency from the 
radio receiving signal, wherein 

the frequency offset estimating step is for determining 
the distortion amounts of the known signals by using the radio 
receiving signal from which the high frequency components have 
been eliminated; and 

the frequency offset eliminating step is for eliminating 
the estimated frequency offset from the radio receiving signal 
whose high frequency components have not been eliminated by the 
filtering step. 

11 . An automatic frequency control method for controlling 
the frequency of a radio receiving signal by eliminating 
frequency offset from the radio receiving signal periodically 
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including a plurality of adjacent known signals, comprising: 

a first distortion amount operating step for determining 
the distortion amounts between the known signals included in 
the radio receiving signal; 

a first frequency offset estimating step for estimating 
the frequency of the direct wave of the radio receiving signal 
based on the distortion amounts of the adjacent known signals 
among the determined distortion amounts; 

a first frequency offset eliminating step for eliminating 
frequency offset corresponding to the estimated frequency of 
the direct wave from the radio receiving signal; 

a second distortion amount operating step for determining 
the distortion amounts of the known signals included in the 
radio receiving signal from which the frequency offset has been 
eliminated; 

a second frequency offset estimating step for estimating 
the center frequency of the Doppler spread of the radio 
receiving signal based on the distortion amounts between the 
periodically included known signal blocks among the determined 
distortion amounts; 

a second frequency offset eliminating step for 
eliminating frequency offset corresponding to the estimated 
center frequency of the Doppler spread from the radio receiving 
signal. 

12 . An automatic frequency control divice for controlling 
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the frequency of a digital baseband signal by eliminating 
frequency offset from the digital baseband signal upon 
receiving the digital baseband signal outputted from an A/D 
converter as an input corresponding to a radio receiving signal 
which periodically includes a plurality of adjacent known 
signals, wherein the frequency of the direct wave of the radio 
receiving signal and the center frequency of the Doppler spread 
of the radio receiving signal are estimated based on the 
distortion amounts of the known signals included in the digital 
baseband signal, and frequency offset is eliminated from the 
digital baseband signal based on both frequencies. 
13. A demodulator, comprising: 

a frequency converter circuit for converting a radio 
receiving signal periodically including a plurality of adjacent 
known signals into an analog baseband signal; 

an A/D converter circuit for converting this analog 
baseband signal into a digital baseband signal; 

a digital signal processing device, which receives the 
digital baseband signal generated by the A/D converter circuit 
as an input, estimates the frequency of the direct wave of the 
radio receiving signal and the center frequency of the Doppler 
spread of the radio receiving signal based on the distortion 
amounts of the known signals included in the inputted digital 
baseband signal, eliminates frequency offset from the digital 
baseband signal based on both frequencies, and eliminates 
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fading distortion from the digital baseband signal from which 
the frequency offset has been eliminated and then demodulates 
the digital baseband signal. 

14. A demodulator according to claim 13, wherein 
the frequency converter circuit has a voltage control 
oscillation part, which oscillates a local oscillation signal 
for converting the radio receiving signal into an analog 
baseband signal, and changes the oscillation frequency of the 
oscillation signal in accordance with an applied voltage; and 
the digital signal processing device eliminates 
frequency offset from the digital baseband signal by applying 
a voltage in accordance with both of the estimated frequencies 
to the voltage control oscillation part. 
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Declaration and Power of Attorney For Patent Application 
Japanese Language Declaration 



As a below named inventor, I hereby declare that: 



My residence, mailing address and citizenship are as 
stated next to my name. 



I believe I am the original, first and sole inventor (if only one 
name is listed below) or an original, first and joint inventor 
(if plural names are listed below) of the subject matter 
which is claimed and for which a patent is sought on the 
invention entitled. 



METHOD AND APPARATUS FOR AUTOMATIC 
FREQUENCY CONTROL AND DEMODULATOR 
( as amended ) 

the specification of which 

□ is attached hereto. 



March ^1, 2000 



as United States Application Number or PCT 
International Application Number 



PCT/JPOO/02123 and was amended on 
(if applicable) 



I hereby state that I have reviewed and understand the 
contents of the above identified specification, including the 
claims, as amended by any amendment referred to above. 



I acknowledge the duty to disclose information which is 
material to patentability as defined in Title 37, Code of 
Federal Regulations, § 1.56. 
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mt. *B^feft^35li119^(a)-(cl):®X}i365^ (b) 5 



Prior Foreign Application(s) 



(Number) (Country) 



1 hereby claim foreign priority under Title 35, United States 
Code, § 119 (a)-(d) or 365(b) of any foreign application(s) 
for patent or inventor's certificate, or § 365(a) of any PCT 
International application which designated at least one 
country other than the United States, listed below and have 
also identified below, by checking the box, any foreign 
application for patent or inventor's certificate, or PCT 
international application having a filing date before that of 
the application on which priority is claimed. 

Priority Claimed 

□ □ 

(Day/MonthA'ear Filed) Yes No 

I hereby claim the benefit under Title 35, United States 
Code, §11 9(e) of any United States provisional 
application(s) listed below. 



(Application No.) (Filing Date) 

fj^m65^ (c) icS^XttfiJ^^iSd^KL^-to ^fc. 



PCT/JPOO/02123 March 31, 2000 

(Application No.) (Filing Date) 



(Application No.) (Filing Date) 

I hereby claim the benefit under Title 35, United States 
Code, § 120 of any United States application(s), or § 365(c) 
of any PCT International application designating the United 
States, listed below and, insofar as the subject matter of 
each of the claims of this application is not disclosed in the 
prior United States or PCT International application in the 
manner provided by the first paragraph of Title 35, United 
States Code, § 112, I acknowledge the duty to disclose 
information which is material to patentability as defined in 
Title 37, Code of Federal Regulations, § 1 .56 which became 
available between the filing date of the prior application and 
the national or PCT International filing date of this 
application. 



(Status: Patented, Pending, Abandoned) 



(Application No.) 



(Filing Date) 



&-5<mmti^xMmxh^ tmzxv^^zh. ^himm 
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(Status: Patented, Pending, Abandoned) 

I hereby declare that all statements made herein of my own 
knowledge are true and that all statements made on 
information and belief are believed to be true; and further 
that these statements were made with the knowledge that 
willful false statements and the like so made are punishable 
by fine or imprisonment, or both, under Section 1001 of Title 
18 of the United States Code and that such willful false 
statements may jeopardize the validity of the application or 
any patent issued thereon. 

POWER OF ATTORNEY: As a named inventor, I hereby 
appoint the following attorney(s) and/or agent(s) to 
prosecute this application and transact all business in the 
Patent and Trademark Office connected therewith: (list 
name and registration number) 
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(703) 413-3000 





Full name of sole or first inventor 
Takashi Asahara. 




Inventor's signature Date 




Residence 

Chiyoda-ku, Tokyo. Japan --J P/ 


mm 


Citizenship 
Japan 




Mailing Address 

c/o MITSUBISHI DENKl KABUSHIKl KAISHA, 2- 
3, Marunouchi 2-chome, Chiyoda-ku, Tokyo 100- 
8310, Japan 






Full name of second joint inventor, If any 

c- lasMiaiu Kojinia. 




Second inventor's signature . Date 




Residence 'f.k.- 
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Japan 
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